Wave propagation in anisotropic elastic materials and curvilinear coordinates using a summation-by-parts finite difference method

نویسندگان

  • N. Anders Petersson
  • Björn Sjögreen
چکیده

We develop a fourth order accurate finite di↵erence method for solving the three-dimensional elastic wave equation in general heterogeneous anisotropic materials on curvilinear grids. The proposed method is an extension of the method for isotropic materials, previously described in the paper by Sjögreen and Petersson [J. Sci. Comput. 52 (2012)]. The proposed method discretizes the anisotropic elastic wave equation in second order formulation, using a node centered finite di↵erence method that satisfies the principle of summation by parts. The summation by parts technique results in a provably stable numerical method that is energy conserving. We also generalize and evaluate the super-grid far-field technique for truncating unbounded domains. Unlike the commonly used perfectly matched layers (PML), the super-grid technique is stable for general anisotropic material, because it is based on a coordinate stretching combined with an artificial dissipation. As a result, the discretization satisfies an energy estimate, proving that the numerical approximation is stable. We demonstrate by numerical experiments that, if the super-grid layers are su ciently wide, the errors due to truncating the domain are of the same order as, or smaller than, the propagation errors from the interior of the domain. Applications of the proposed method are demonstrated by three-dimensional simulations of anisotropic wave propagation in crystals.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

An efficient finite difference time domain algorithm for band structure calculations of Phononic crystal

In this paper, a new algorithm for studying elastic wave propagation in the phononic crystals is presented. At first, the displacement-based forms of elastic wave equations are derived and then the forms are discretized using finite difference method. So the new algorithm is called the displacement-based finite difference time domain (DBFDTD). Three numerical examples are computed with this met...

متن کامل

MPI- and CUDA- implementations of modal finite difference method for P-SV wave propagation modeling

Among different discretization approaches, Finite Difference Method (FDM) is widely used for acoustic and elastic full-wave form modeling. An inevitable deficit of the technique, however, is its sever requirement to computational resources. A promising solution is parallelization, where the problem is broken into several segments, and the calculations are distributed over different processors. ...

متن کامل

Wave Propagation in Sandwich Panel with Auxetic Core

Auxetic cellular solids in the forms of honeycombs and foams have great potential in a diverse range of applications, including as core material in curved sandwich panel composite components, radome applications, directional pass band filters, adaptive and deployable structures, filters and sieves, seat cushion material, energy absorption components, viscoelastic damping materials and fastening...

متن کامل

Simulation of Wave Propagation Along Fluid-Filled Cracks Using High-Order Summation-by-Parts Operators and Implicit-Explicit Time Stepping

We present an efficient, implicit-explicit numerical method for wave propagation in solids containing fluid-filled cracks, motivated by applications in geophysical imaging of fractured oil/gas reservoirs and aquifers, volcanology, and mechanical engineering. We couple the elastic wave equation in the solid to an approximation of the linearized, compressible Navier–Stokes equations in curved and...

متن کامل

Significant Error Propagation in the Finite Difference Solution of Non-Linear Magnetostatic Problems Utilizing Boundary Condition of the Third Kind

This paper poses two magnetostatic problems in cylindrical coordinates with different permeabilities for each region. In the first problem the boundary condition of the second kind is used while in the second one, the boundary condition of the third kind is utilized. These problems are solved using the finite element and finite difference methods. In second problem, the results of the finite di...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • J. Comput. Physics

دوره 299  شماره 

صفحات  -

تاریخ انتشار 2015